Discrepancy between hadron matter and quark-gluon matter in 

net charge transfer fluctuation 

Dai-Mei Zhou 1 , Xiao-Mei Li 2 , Bao-Guo Dong 2 , and Ben-Hao Sa 3 ' 2 - 1 ' 4 * 
1 Institute of Particle Physics, Huazhong Normal University, Wuhan, 430079 China 

2 China Institute of Atomic Energy, 
P. 0. Box 275 (18), Beijing, 102413 China 
3 CCAST (World Lab.), P. O. Box 8730 Beijing, 100080 China 
4 Institute of Theoretical Physics, Academy Sciences, Beijing, 100080 China 

Abstract 

A parton and hadron cascade model, PACIAE, is employed to investigate the net charge 
transfer fluctuation within \n\=l in Au + Au collisions at y / s]v r /v=200 GeV. It is turned out that 
the observable of net charge transfer fluctuation, k , in hadronic final state (HM) is nearly a factor 
of 3 to 5 larger than that in initial partonic state (QGM). However, only twenty percent of the 
net charge transfer fluctuation in the QGM can survive the hadronization. 
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Recently the net charge transfer fluctuation has been proposed in cLS cL SI enal of 

Quark-Gluon-Plasma (QGP) phase transition expected to be existing in the relativistic 
nucleus-nucleus collisions. Following 3j the observable 

k( V ) = D u ( V )/(dN ch /d V ) (1) 

is employed to describe the net charge transfer fluctuation. In the above equation the net 
charge transfer deviation, D u , reads 

DM =< u(7]) 2 >-< u( V ) > 2 , (2) 

where the net charge transfer, u, is defined by 

u( V ) = [Q F ( V )-Q B ( V )]/2, (3) 

where the Qf{v) {Qb{v)) is referred to the net charge in forward (backward) region of rj 
and the N c h stands for the charge multiplicity accounted according to the charge of particle. 
The k is argued to be a measure of the local unlike-sign charge correlation length P and 
the charge correlation length in QGP phase (in quark-gluon matter, QGM) is expected to 
be much smaller than the one in hadronic matter (HM) because the charge unit is 1/3 and 
1 in QGM and HM QQEJ, respectively. 

In |2j a neutral cluster model was used first and the hadronic transport models of HIJING 
RQMD [^|, and UrQMD [l| were employed then to study the net charge transfer fluctu- 
ation. The results from above hadronic transport models could be summarized as follows: 
1. The discrepancy among them is not obvious from each other. 2. The k{t) = 1) calculated 
in interval of |?y| < 1 is equivalent to the net charge fluctuation at r\ — 1 and is close to 
the STAR datum fljl of ~ 0.27 ± 0.02 (cited from Q directly) in Au + Au collisions at 
y/sNN=130 GeV. 3. The k(t]) does not strongly depend on the centrality. 

A parton and hadron cascade model, PACIAE, is employed in this letter investigating 
the net charge transfer fluctuation, k, within \q\ < 1 both in the early partonic stage (QGM) 
and in the hadronic final state (HM) in Au + Au collisions at v /saT/v=200 GeV. As expected 
the later results are quite close to the results in HIJING, RQMD, and UrQMD. However 
the former results are smaller than the later one by a factor of 3 to 5. Unfortunately, the k 
in QGM seems to be hard to survive the hadronization. 

As the simplified version of PACIAE model has been published in and a nice bit of 
detailed description has been given in [l^], here we just give a brief introduction for PACIAE 



model. In the PACIAE model a nucleus-nucleus collision is decomposed into nucleon-nucleon 
collisions. The nucleons in a nucleus is distributed randomly according to Wood-Saxon 
distribution. A nucleon-nucleon collision is described by the Lowest-Leading-Order (LLO) 
pQCD parton-parton hard interactions with parton distribution function in a nucleon and by 
the soft interactions considered empirically, that is so called "multiple mini-jet production" 
in HIJING model j?|. However, in PACIAE model that is performed by PYTHIA model 
with string fragmentation switched-off. Therefore, the consequence of nucleus-nucleus 
collision is a configuration of a (q), diquark (anti-diquark), and g, besides the spectator 
nucleons and beam remnants [13j . The diquark (anti- diquark) is forced to split into qq (qq) 
randomly. 

So far we have introduced the partonic initialization of nucleus-nucleus collision in PA- 
CIAE model, what follows is then parton evolution (scattering). To the end, the 2— >2 LLO 



pQCD differential cross section 



14( is used. Of course, that must be regularized first by 



introducing the color screen mass. The total cross section of parton i bombarding with j 
could then be calculated via a integral over the squared momentum transfer in a subprocess 
ij — > kl and a summation over partons k and I. With above differential and total cross 
sections the parton scattering can be simulated by Monte Carlo method. As of now, only 
2— >2 processes are involved, among them there are six elastic and three inelastic processes 

0. 

As for the hadronization we first assume that the partons begin to hadronize when the 
interactions among them have been ceased ( freeze-out). They could hadronize by either 
Mentation .node, QQ or coatee m „de, QQ. What the .mentation mod e.s 
included here are the Field- Feynman model, i. e. Independent Fragmentation (IF) model 



and Lund string fragmentation model |l6j. However, the program built in is employed 
for the implementation of fragmentation model. On the contrary, we do write a program for 
coalescence model ourselves. 

The hadron evolution (hadronic rescattering) is modeled as usual two body collisions and 
is copied directly from a hadron and string cascade model LUCIAE [3]. There is no need 
to say more about the hadronic rescattering referring to ljj if necessary. 
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FIG. 1: Observable of net charge transfer fluctuation, k, as a function of r] within \rj\ < 1 in 0-5% 
most central Au + Au collisions at ^/saw=200 GeV. 

Since we are not aimed to reproduce the experimental date but to study the physics of net 
charge transfer fluctuation, we do not adjust model parameters at all. In the calculations the 
IF model is adopted for hadronization and the net charge transfer fluctuation is counted 
in the interval of \r)\ <1. The simulated results by default PACIAE model are indicated with 
"HM w/ QGM" (HM with QGP assumption), since the hadronic final state is evolved from 
partonic initial state. If the simulation is ended up at the stage of partonic scattering and 
the net charge transfer fluctuation is counted over partons only, the results will be referred 
to as "QGM" . In that calculation it is assumed that the gluon does not contribute to the net 
charge but it does contribution to charge multiplicity by 2/3 as assumed in 0, ^ the 
simulation is ended up at the stage of partonic scattering and both the partons and beam 
remnants (hadrons) are counted in net charge transfer fluctuation the corresponding results 
are then symboled as "QGM w/ remnant". It should be mentioned here that the spectator 
nucleons do not affect the net charge transfer fluctuation in \rj\ <1. A calculation where 
the string fragmentation in PYTHIA is switched-on and followed directly by the hadronic 
rescattering is referred to as "HM", since in this simulation only hadronic transport is taken 



221. 



into account, like in HIJING, RQMD, UrQMD, and JPCIAE 

Fig. Ogives the observable of net charge transfer fluctuation, k, as a function of pseudo- 
rapidity, rj, in the simulations of "HM", "HM w/ QGM", and "QGM" (solid circles, squares, 
and triangles, respectively) in 0-5 % most central Au + Au collisions at v /J/v r /v=200 GeV. 
One sees in this figure that the k of "HM w/ QGM" reproduces nicely the STAR datum and 
the k of "HM" is a bit larger than the STAR datum at rj—1. The trend of k varying with 
77, both in "HM" and "HM w/ QGM", is similar to the ones in HIJING, RQMD, UrQMD 
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(cf. Fig. 5 in Q]). On the contrary, the k of "QGM" keeps nearly constant, like the charge 
fluctuation as function of rapidity interval in QGM in thermal model "\ and in transport 
model 2lj]. It is interesting to see that the k in "HM" is larger than k in "QGM" by a factor 
of 3 to 5 from upper rj to the lower rj. The discrepancy between k in "HM" and in "HM 
w/ QGM " amounts ~ 20% in the average, that means that the probability of net charge 
transfer fluctuation in QGM surviving hadronization can be estimated to be ~ 20% either. 
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FIG. 2: k as a function of ry within \rj\ < 1 in a range centralities of the Au + Au collisions at 
^/i^7=200 GeV. 

In Fig. |21 the centrality dependence of k(t/) is given both for "HM" and "QGM" (solid 
and open symbols, respectively) for 0-5, 30-40, and 70-80% central Au + Au collisions at 
v /s^7=200 GeV. Both of the k(tj) in "HM" and in "QGM" do not strongly depend on the 
centrality which is consistent with the results in hadronic transport models HIJING, RQMD, 
and UrQMD (cf. Fig. 5 in J)- 

We compare the k(t]) in "QGM" to the one in "QGM w/ remnant" in Fig. 01 for 0-5% 
most central Au + Au collisions at y / J/v r /v=200 GeV. One sees in this figure that the influence 
of beam remnants upon the k(t]) in "QGM" amounts ~35% in average. This influence does 
not change the status of big difference between k in "HM" and in "QGM" shown in Fig. 
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FIG. 3: k as a function of rj within |?7| < 1 in 0-5% most central Au + Au collisions at y / s/v r iv=200 
GeV. 

In summary, a parton and haclron cascade model, PACIAE, is employed investigating the 
net charge transfer fluctuation within |?7|=1 both in partonic initial state and in hadronic 
final state for a range centralities of Au + Au collisions at v /J/v r /v=200 GeV. In the hadronic 
final state (k in "HM") the observable of net charge transfer fluctuation, k, turns out to be 
nearly a factor of 3 to 5 larger than the k in partonic initial state (k in "QGM"). However, 
the k in "QGM" is hard to survive the hadronization, the survival probability amounts ~ 
twenty percent. 
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